Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN I^TeX style file v2.2) 



Irradiation of accretion discs in active galactic nuclei due 
to warm absorber 



Z. Loska\ B. Czerny^, R. Szczerba^ 

^ N. Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 
^N. Copernicus Astronomical Center, Rabianska 8, 87-100 Torun, Poland 



2 February 2008 



o 
o 

(N 

CD 



(N 

o 
o 



ABSTRACT 

The presence of the warm absorber of considerable optical depth is seen in many AGN. 
We show that this medium may affect significantly the optical/UV spectrum of an 
AGN by backscattering a fraction of the total radiation flux towards the disc surface. 
We consider in detail the case when the disk extends down to a marginally stable 
orbit, all the emission comes from the disk surface and the scattering medium forms 
a cone around the symmetry axis. Disc irradiation results in much flatter optical/UV 
continuum than predicted by standard disc models. The effect depends both on the 
total optical depth of the warm absorber and on the specific density distribution of this 
medium so the analysis of the optical/UV continuum allows to obtain constraints for 
the warm absorber complementary to those obtained from the soft X-ray data analysis. 
We give results for two exemplary sources - RE J1034-I-396 and PG1211-I-143, and for 
the bluest composite quasar spectrum of Richards et al. obtained from SDSS. 

Key words: Radiative transfer, Accretion discs, Galaxies:active, Galaxies:Seyfert, 
X-rays:galaxies 
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1 INTRODUCTION 

Optical/UV spectra of bright active galactic nuclei (AGN) 
are usually interpreted as originating from an accretion disc 
surrounding the central massive black hole (Lynden-Bell 
1969, Shields 1978, Malkan & Sargent 1982; for a review, 
see e.g. Koratkar & Blaes 1999, Czerny 2003). However, both 
the shapes of the radiation spectra and the variability stud- 
ies (e.g. Wanders et al. 1997, Collier et al. 1998 for NGC 
7469; see Ulrich, Maraschi & Urry 1997 for a review) sug- 
gest that significant fraction of radiation emitted at a given 
radius is due to the reprocessing of the radiation dissipated 
elsewhere. 

Most of the radiation is produced in the deep potential 
well, close to the black hole. In quasars and Seyfert 1 galaxies 
part of this radiation reaches an observer directly. However, 
part of the radiation may be recaptured by the outer disc 
region. 

The irradiation of the disc may be direct or indirect. 
In the first case the disc must be exposed to the radiation 
fiux generated in the inner region. In the second case there 
may be a scattering medium which redirects a part of the 
radiation back towards the disc. 

The direct irradiation of AGN accretion discs was dis- 
cussed by several authors (e.g. Krolik et al. 1991, Cour- 
voisier & Clavel 1991, CoUin-Souffrin 1991, Rokaki et al. 
1993, Loska & Czerny 1997, Kurpiewski et al. 1997, Sincell 



& Krolik 1997, Collin & Hure 1999, Nayakshin et al. 2000, 
Ballantyne et al. 2001, Nayakshin & Kazanas 2002, Soria & 
Puchnarewicz 2002, Rozanska et al. 2002, Chiang & Blaes 
2001, 2003, Chiang 2002). Most of these papers were actu- 
ally devoted to the irradiation of the inner disc by the hard 
X-rays generated in a lamp-post or disc corona. Much of the 
research had been done in the context of X-ray binaries (e.g. 
Raymond 1993, Dubus et al. 2001, Jimenez-Garate et al. 
2002) or cataclysmic variables (Meyer & Meyer-Hofmeister 
1982, Smak 1989) and the methods are directly applicable 
to AGN discs but the results may not apply quantitatively 
since the disc shape depends on the mass of the central ob- 
ject, as already seen from the classical paper of Shakura & 
Sunyaev (1973). 

The indirect irradiation was also considered in several 
papers (Ostriker et al. 1991, Murray et al. 1994, Kurpiewski 
et al. 1997; also Esin et al. 1997 in the case of galactic 
sources). These authors analyzed the scattering of the ra- 
diation in the accretion disc corona. 

In the present paper we consider the scattering of the 
disc radiation by the warm absorber. The warm absorber 
medium, seen in many AGN (for a review, see e.g. Kriss 
2004), is a highly ionized plasma so it predominantly scat- 
ters the passing radiation, with absorption important only in 
the soft X-ray band. The optical depth of the warm absorber 
for scattering is difficult to measure but estimated values 
can be as high as 0.33 or more for PG1211-I-143 (Pounds 
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et al. 2003a), 0.26 or more for PG PG 0844+349 (Pounds 
et al. 2003b), 0.20 or more for PG 1402+261 (Reeves, Por- 
quet & Turner 2004). Arav et al. (2003) show by analyzing 
OVI UV absorption and OVI X-ray absorption in NGC 5548 
that simple modelling gives discrepant results and in general 
heavily underestimates the amount of outflowing material. 

Therefore, a significant fraction of the radiation flux 
emitted by the central region can be in such way redi- 
rected towards the outer parts of the accretion disc and 
thermalized there. Therefore, the presence of highly ionized 
warm absorber of considerable optical depth must modify 
the IR/optical/UV disc spectra. We can apply this effect ei- 
ther to explain the departure of the observed spectrum from 
the simple prediction of the accretion disc model or, if no 
such clear departure is seen, as an independent way to put 
limits to the optical depth and the covering factor of the 
warm absorber. 

The effect should be present in all AGN, but it is po- 
tentially most important when the direct disc irradiation by 
an X-ray source is relatively weak. We assume the accre- 
tion flow proceeds through a standard disc which extends 
down to the marginally stable orbit. We neglect any hard 
X-ray emission. We calculate the optical/UV spectra of the 
accretion discs taking into account the scattering of the disc 
flux by the warm absorber as well as the effect of the direct 
self-irradiation of the disc. 

As an application, we show the spectra models for two 
NLSl: RE J1034+396 and PG 1211+143, and for the com- 
posite 1 quasar spectrum of Richards et al. (2003). 



2 METHOD 

We consider the self-irradiation of the disc both due to di- 
rect self-illumination and due to radiation scattering by the 
warm absorber. 

In the case of an extended scattering medium, there is 
practically no dependence of the result on the disc shape. 
However, the direct irradiation is very sensitive to the 
adopted disc shape: the disc must flaring in order to have its 
surface exposed to the central source. We must specify the 
disc model appropriate for an AGN in order to see which 
disc parts actually are irradiated. 

Therefore, in Sect. O we describe the disc model, in 
Sect. 12.21 we show the formulae for the direct irradiation, 
in Sect. 17!^ we describe the warm absorber model and the 
method of computing the incident flux at a given disc radius, 
and in Sect. 12.41 we describe the method of calculating disc 
spectrum. 



2.1 Disc shape 

We use the stationary Keplerian disc model of Rozaiiska et 
al. (1999). The disc thickness is calculated by solving the 
equations of the disc vertical structure. Viscosity is described 
assuming the standard a approach of Shakura & Sunyaev 
(1973), i.e. we assume that the viscous torque is propor- 
tional to the total (gas + radiation) pressure. Recent stud- 
ies of magneto-rotational instability generally supports the 
validity of such scalling (e.g. Turner 2004). Estimates based 
on the MHD simulations of the MRI instability indicate the 
average value of the effective viscosity parameter a about 



0.02 (Winters et al. 2003), and similar values of a between 
0.01 and 0.03 were derived from the analysis of the AGN 
variability (Starling et al. 2004). 

The disc vertical structure (as a function of the distance 
2 from the disc equatorial plane) is calculated by solving the 
equations of viscous energy dissipation, hydrostatic equilib- 
rium, and energy transfer: 
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Here is the Keplerian angular velocity, a and c are physical 
constants. Fi is the energy flux transported locally in the 
direction perpendicular to the equatorial plane and carried 
either by radiation or by convection: 



Fl = Frad 

Fl — -^rad + Fcon 



Vrad < Vad 
Vrad > Vad 



(4) 



These equations, supplemented by the equation of state, de- 
scribe the temperature, density and pressure proflles, T{z), 
p{z) and P{z). 

The frequency-averaged opacity k is taken to be the 
Rosseland mean and includes the electron scattering, free- 
free and bound-free transitions. The opacity tables are from 
Alexander, Johnson & Rypma (1983) and Seaton et al. 
(1994). The presence of dust and molecules is included in 
the opacity description. The details of the model were dis- 
cussed in Pojmariski 1986 and Rozaiiska et al. (1999). 

Disc thickness is determined iteratively, from the con- 
dition that the energy flux dissipated inside the disc is equal 
to the energy flux given by assumed accretion rate through 
the standard formula 



F{H,r) 
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This means that we use the Newtonian approximation for 
a non-rotating black hole. Consequently, we use the dimen- 
sionless accretion rate in Eddington units, m = M/Mfidd, 
with the efficiency of 1/12: 



Me 



2.66 



M 



108 Mf 



-[Me/yr]. 



(6) 



The effect of irradiation is included only through the 
additional term present in Eq. |K| (see e.g. Tuchman et al. 
1990), where Fine is either Fdirect or Fsc, or both, since we 
concentrate on the overall spectral shape, particularly in the 
optical/UV band. This term describes the absorbed fraction 
of the external incident flux. The method of determination 
of this flux is described in the Sections 12.21 and 12.31 The 
albedo. A, is assumed to be independent from the photon 
energy and its energy-independent mean value is assumed 
to be equal to 0.2, appropriate for cold disc surface (see e.g. 
Haardt & Maraschi 1991). 



2.2 Direct irradiation 

Computing the direct irradiation we follow the approach of 
Fukue (1992). 
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Figure 1. The schematic picture of the self-irradiation of the 
accretion disc. 



We calculate the irradiation flux from the given distri- 
bution of the emitted flux, Fem{r) and the disc shape, H{r). 
The incident flux at a given disc radius, Vo, has the form of 
the integral over the disc surface 

4 

^direct (^o) 

>. dr ' 



Rout F(H,r) cosicosio J 1 + {^Yrdrd(t) 



(7) 

/O J Ri, 

with an additional condition that the contribution to the 
integral is included only if 

cosi>0, cosio>0. (8) 

The incident flux Fine ~ Fdirect is calculated at the ra- 
dius To, r is the current radius of integration over the disc 
surface. Ho is the disc thickness at Vo while H is the disc 
thickness at the radius r, io and i are the corresponding 
angles between the normal to the disc surface and the line 
joining two disc surface points, as marked schematically in 
Fig. Q The distance Ds between the irradiating and irradi- 
ated surface element is given by 



Dt 



T +ro + {H-Ho) -2rroC0S( 



the angle i is determined by 

cos I — — — , 

and the angle io is determined by 

^{ro -r cos rt>)+H - Ho 

COSlo = — = , 



(9) 



(10) 



(11) 



Radiation flux from the disc is determined by the Eq.|S| 
Since the effect of irradiation depends on the shape of the 
disc, and both the disc shape and the local emissivity depend 
on irradiation, the result is obtained through consecutive 
iterations, starting from unilluminated disc. Usually 2-3 
iterations are enough to provide the result with the accuracy 
better than 1% . 



2.3 Scattering by the warm absorber 

The idea of redirecting some of the disc radiation flux emit- 
ted in the central region towards the outer part through 
scattering by the warm absorber was suggested by Czerny 



et al. (2003). It could be seen from Fig. 3 of this paper that 
the warm absorber absorbs efficiently only in the soft X-ray 
band while for lower frequency photons the opacity is equal 
to that of electron scattering and this effect dominates the 
overall extinction in the optical/UV band. The key param- 
eter is the ionization parameter ^ defined as 



L 



(12) 



where L is the luminosity of the central source, n is the 
number density of the shell and R is the inner radius of the 
shell. 

We checked the effect using the code CLOUDY, version 
94.0, for a spherical shell geometry. The reflected contin- 
uum is redirected towards the disk and the transmitted con- 
tinuum is seen by an observer. For ^ — 10"* the warm ab- 
sorber imprints significant features in the X-ray band into 
the transmitted spectrum which allow to determine the hy- 
drogen column of the absorbing medium. Narrow absorption 
lines form also in the UV band. At the same time the overall 
reflected spectrum still is not much different from the intrin- 
sic continuum. In the case of almost fully ionized absorber 
(^ — 10^) electron scattering dominates so strongly that 
actually almost no spectral features are seen. It supports 
even better our later approach based on the negligence of 
absorption processes. At the same time, it illustrates poten- 
tial problems in detecting warm absorbers. An extremely 
strongly ionized medium with large optical depth may well 
exist in an object and still escape detection by analysis of 
the soft X-ray data. 

The presence of the dust within the warm absorber com- 
plicates the issue. As an example, we considered a case with 
C = 10", but we allowed for an additional fraction of dust, 
at amount of 1% of the standard gas to dust ratio for the 
interstellar medium in our Galaxy, and of the same com- 
position. Electron scattering is still quite effective, and the 
dust also partially scatters the incident photons. The net 
effect is strong reddening in the transmitted continuum and 
relatively minor suppression of the UV part in the reflected 
spectrum. Details of those effects depend significantly on the 
dust composition. The possible presence of the dust in the 
warm absorber is discussed in a number of papers (e.g. Ko- 
mossa & Breitschwerdt 2000, Mason et al. 2003, Ballantyne, 
Weingartner & Murray 2003) . 

In the present paper we concentrate on objects which 
are not heavily reddened in UV and we consider the warm 
absorber medium as purely scattering fully ionized plasma. 

The origin, the geometry and the physical state of the 
warm absorber are still under discussion. As for the source 
of material, ideas include accretion disk winds/Broad Line 
Region clouds (Elvis 2000, Kriss et al. 2003, Blustin et al. 
2003), inner edge of the dusty /molecular torus (Krolik & 
Kriss 2001) and radiation driven clouds constituting the 
Narrow Line Region (Crenshaw et al. 2000). King & Pounds 
(2003) go as far as suggesting that actually the outflow from 
the direct vicinity of an accreting black hole may be strong 
enough to form an optically thick shell that Comptonize the 
emission from the directly unseen inner source. 

Narrow absorption lines in many spectra of Seyfert 1 
galaxies indicated a clumpy medium with numerous colder 
clumps embedded in an fully ionized hotter medium (Cren- 
shaw et al. 1999; Kriss 2002; Krolik & Kriss 2001, Blustin et 
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Figure 2. The schematic picture of scattering of the disc emission 
by the warm absorber. 



al. 2003) while 280 ksec XMM-Newton spectroscopy of NGC 
3783 (Behar et al. 2003) suggested mostly uniform medium, 
at distances up to 2.8 pc from the nucleus. 

The optical depth of the warm absorber is poorly con- 
strained since the directly measured column density of spe- 
cific ions must be supplemented with modelling of the ion- 
ization state of the plasma in order to obtain the hydrogen 
column. 

If cooler clumps are embedded in a fully ionized medium 
the contribution of this medium is particularly difficult to 
estimate while it may significantly enhance the total optical 
depth of the medium for electron scattering. 

Detection of the variability in the OVIII hydrogen col- 
umn of MCG-6-30-15 (Otani et al. 1996) allowed to estimate 
the distance of the warm absorber as smaller than 10^^ cm 
but the lack of variability in the OVII edge suggested that 
the warm absorber is rather spatially extended, with highly 
ionized fiow closer in and less ionized flow further out, as 
visualized by Blustin et al. (2003). Models based on evapo- 
ration of the inner torus suggest distances of a fraction of a 
parsec while models with disk outflow suggest distances of 
a few Schwarzschild radii for fully ionized flow. 

Radial velocities of warm absorber features of order of 
a few hundreds km/s determined in a number of sources, if 
interpreted as roughly Keplerian speed, may indicate rather 
large distances and the connection with dusty torus (e.g. 
Ashton et al. 2004). Large velocities claimed to be seen 
in PG1211-H143 (Pounds et al. 2003a) and PG0844-h349 
(Pounds et al. 2003b) suggest much smaller values, although 
identification of the absorption lines and consequently, de- 
termination of the outfiow velocity, may be questioned 
(Kaspi 2004). 

We assume that the warm absorber has an axially- 
symmetric conical shape. It extends from the inner radius, 
Rmin, to the outer radius, Rmax, with the opening angle of 



the cone given by Omax- The schematic picture is shown in 
Fig.ia 

The density of the warm absorber depends on the radial 
distance from the center, R in the form 



n{R)^n^{-^)-^ 



(13) 



where the density at the inner edge of the warm absorber is 
conveniently expressed as 

(14) 



-l^dR 



since the total optical depth of the warm absorber, rtot, is 
relatively easily estimated from observations. We can also 
include in our scheme the dependence of the density distri- 
bution on 9. 

Since the warm absorber is an extended medium, we ne- 
glect the extension of the inner parts of the disc where most 
of the energy is generated and we assume that the whole 
radiation comes from a point-like source. Its luminosity, L, 
is determined by the accretion rate 



L = mLEdd- 



(15) 



We take into account the finite optical depth of the 
warm absorber so the local emissivity of this medium due 
to the electron scattering is given by 

j{R,e) = aT{R) '2L^ ^ (^g) 

47r7?2 exp[— aT{s)ds\ 

where the opacity ar, is the Thompson opacity [axiR) = 
(TTn(R)). 

This extended medium returns the fraction of the scat- 
tered radiation towards the disc. The resulting illumination 
is calculated as 



Fsc{ro) 



^ r2Ti rR 
47r Jo Jr^ 
where 



j{R, 6)R^ sin 9 cos io 

-'-^11} 



dedRdfj), (17) 



cos to 



and 



dHg 
dr 



{vo — R sin 9 COS tp) + R cos 9 — Ho 



Di = rl + R^ + Ho - 2R{ro sin 61 cos <^ -I- Ho cos 9). 



(18) 



(19) 



2.4 Accretion disc spectrum 

Accretion discs in high Eddington rate AGN are opti- 
cally thick and they thermalize the absorbed radiation effi- 
ciently. The models of the emitted local spectrum, however, 
show some departure from the local black body approxi- 
mation due to the role of the Compton scattering in the 
disc outer layers (e.g. Shimura & Takahara 1995, Madej & 
Rozanska 2004). This departure is usually well described 
in terms of the colour temperature to effective tempera- 
ture ratio. In various models this ratio varies from 1.7 - 2.0 
(Shimura & Takahara 1995, Ross & Fabian 1996) to 1.8 - 2.6 
(Merloni, Fabian & Ross 2000). Madej & Rozanska (2004) 
find Tcoi/Teff < 2 in their set of neutron star atmosphere 
models. Recent observational studies of galactic sources in 
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Figure 3. The dependence of the disc thickness on the disc radius 
without the effect of self-irradiation (dashed hne), with the self- 
irradiation (continuous line) and with the back-scattering flux by 
the warm absorber (dotted line). Disc parameters: M = 10^ Mq, 
m = 0.5, a = 0.02, warm absorber parameters: rtot = 0.6, /3 = 0, 

Rmin = WOORschw J Rmax = ^'^OORschw t 9max = 7I"/3. 



their high states (Gierliriski & Done 2004a) suggest that the 
value of ~ 1.8 indeed well represents the effect. We adopt 
this value in our description of the inner part of the disc 
where electron scattering dominates the opacity. In the out- 
ermost region the disc is not so strongly ionized and the 
emission is expected to be close to a black body emission. 
Therefore, we describe the local spectrum as 



cost 
47r_D2 



■RB,{Tcoi)2Tvrdr, 



(20) 



where the effective temperature is calculated from the total 
(dissipated -I- thermalized) flux 

n4 



eff 



F{H,r) 



(21) 



and to ensure the continuity between the inner and the outer 
disc region, we assume 

Tcoi/T,ff = 2.228 - 0.041 (r/i?sch») if r < 30iisch» (22) 
and Tcoi ~ T^ff at larger radii. 

3 RESULTS 
3.1 Disc shape 

An exemplary shape of AGN accretion disc is shown in 
Fig. 13 It is much more complex than the simple prediction 
of Shakura & Sunyaev (1973) for radiation- pressure domi- 
nated region. The departure is not important from the point 
of view of irradiation by the extended warm absorber but it 
is important for self-irradiation. Even in the innermost ra- 
diation pressure dominated region the disc thickness is not 



constant but rises slowly. Further out, it rises faster as the 
radiation pressure drops. Some wiggling is also visible - this 
is due to the complex opacities used for the description of 
the vertical structure. The opacity of the outermost weakly 
ionized disc region is dominated by molecules and absorp- 
tion by dust, so the disc is never optically thin. 

In particular, the drop in the disc height at 
SOOORschm happens in the region of ionization instability 
(Lin & Shields 1986) . Here we neglect the issue of the time- 
dependent disc behaviour. It is most probably justified for 
high accretion rate objects since no evaporation of an inner 
disc is expected and ionization instability reduces to small 
flickering for a constant viscosity law (Siemiginowska, Cz- 
erny & Kostyunin 1996, Janiuk et al. 2004). If the irradia- 
tion by the warm absorber is included, this region is shifted 
outward (see Fig. |3 dotted line). The disc thickness is in- 
creased by a factor of 3 at distances ~ lOOORschw from the 
black hole since irradiation was most efficient there for an 
adopted model of the warm absorber. Self-irradiation (see 
Fig.|5Jl has a very small impact on the disc thickness, visible 
only at very large distances, above 10''' Rschw 

If accretion rate is significantly lower than used in 
Fig- El i-6- ^ ~ 0-05 instead of m = 0.5 then the disc thick- 
ness in the innermost (radiation pressure dominated) part is 
roughly by a factor 10 smaller, the wiggle in the disc shape 
(i.e. partial ionization instability region) moves in roughly 
by a factor of 3, and the disc thickness in the outer region is 
almost unchanged. The effect of the self-irradiation and irra- 
diation due to the scattering by the warm absorber modify 
the disc thickness in roughly the same way as before. 

The change in the properties of the warm absorber is 
reflected in the disc thickness. Lower/higher total opacity 
leads to smaller /larger increase in the disc height at dis- 
tances where the irradiation is important. 

The density profile of the warm absorber specifies at 
which distance the irradiation is most effective {Rmin, Rmax 
or a broad range between the two). However, performing 
subsequent iterations between irradiation and the disc thick- 
ness affect the final irradiation flux rather weakly. 



3.2 Radial dependence of the irradiating flux 

3.2.1 Direct irradiation 

An example of the incident radiation flux in the absence 
of the warm absorber is shown in Fig. 01 (continuous line). 
Model parameters are the same as used in Fig. |5] 

The flux is large at small radii, but it is nevertheless 
small in comparison with the locally dissipated flux in this 
region. Self-irradiation does not create energy, it only effec- 
tively modifies slightly the angular dependence of the radia- 
tion leaving the central region. However, as it is well known, 
irradiation becomes important at large distances (e.g. Hoshi 
& Inoue 1988). The sudden jump in the fiux is again a con- 
sequence of the rise of the disc thickness at the end of the 
partial ionization zone (see Fig. |3J . 

Sometimes, for simplicity, the irradiation is calculated 
assuming the point-like source instead of an extended inner 
disc. Irradiation flux is in this approximation given by (e.g. 
King 1998, King & Ritter 1998) 



F — 



(23) 
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Figure 4. The irradiation flux (left panel) and the ratio of the irradiation flux to the locally dissipated flux for the same model as 
in Fig. 121 Continuous line - self-irradiation; dotted line - back-scattering by the warm absorber. We also show the back-scattered flux 
(dashed line) obtained for the warm absorber extending from Rmin = ^OORschw to Rmax = 5000Rgchwj a-nd with the same optical 
depth Ttot = 0.6 as before. Thick horizontal line marks the region where irradiation dominates. 



0.1 - > 




1 2 3 4 5 

i°g(R/Rschw) 



Figure 5. The dependence of the 7 factor on the disc radius. Con- 
tinuous line: illumination by extended disc. Dashed line: point-like 
source approximation given by EQ. I24l for n=l. Long dashed line: 
point-like source approximation given by Eg. 1241 for n=2. Model 
parameters as in Fig. |3 Negative 7 means that a point-like source 
located at the center would not irradiate the disc surface. 

where 

7 = (V(^-1). (24) 
r amr 

The form of Eg. 1231 is general so it is illustrative to ex- 
press the distribution of the irradiation flux determined from 
our model also as a radial dependence of the dimensionless 



quantity gamma. The value n = 1 is expected to be better 
for neutron star systems while n = 2 for black hole systems 
(Shakura & Sunyaev 1973, Fukue 1992, King et al. 1997). 

In Fig. |5] we show both the net 7 factor which results 
directly from our computations and the 7 factor from Eg. 1241 
obtained from the disc thickness distribution in our model. 

We see that the irradiation flux computed properly from 
the model is always positive, as expected, since the disc is 
irradiated at least by the nearest region. However, in these 
regions which are not well exposed to the innermost parts 
the irradiation is very small. Effective 7 factor for m = 0.5 
is only ~ 2.5 x 10~^ in the intermediate disc region, between 
100 and 10000 i?schw. It becomes larger (~ 0.0005) at the 
distances 10'' — 10^ ^?Schw, as well as in the innermost region, 
approaching ~ 0.04 at a few i?schw 

Analytical formula (Eq. 1241 dashed line in Fig. |KJ does 
not provide a good approximation of the numerical results. 
Adopting n = 2 in this formula does not lead to better 
approximation - in this case the flux is generally underesti- 
mated, apart from the innermost region of a few Rschv 

Self-irradiation of optically thick, geometrically thin 
disc is generally inefflcient in AGN, as argued already by 
Tuchman et al. (1990). 



3.2.2 Irradiation due to warm absorber 

An example of the incident radiation flux due to warm ab- 
sorber is shown in Fig. 31 (dotted line) . Here we assumed 
the following warm absorber parameters: Rmin = 1000, 
Rm,ax = 1200, Qmax = 7r/3, /3 = 0, Ttot = 0.6. lu this case the 
irradiating flux is almost constant up to r « 1000 and falls 
down as r~^ further out. Therefore, the warm absorber effec- 
tively acts as a single point-like source which captures and 
re-emits the fraction of the disc luminosity approximately 
given by 2r cos and the integration over the volume is not 
essential. 

However, if we assume much larger extension of the 
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warm absorber, the effect will depend essentially on the den- 
sity profile. 

If the power law index describing the radial density pro- 
file, /3, is close to zero the irradiation fiux between Rmin and 
Rmax decreases roughly as r~^, and for /9 ~ 2 decreases as 
. Below Rmin the irradiation flux is constant and above 
Rmax decreases as . 

Generally, the irradiation by the warm absorber is neg- 
ligible in the innermost part of the disc. Self-irradiation is 
relatively more important but not essential, in comparison 
with the radiation fiux dissipated inside the accretion disc. 
At r « Rmin, the irradiation by the warm absorber be- 
comes dominant over self-irradiation and at a few hundred 
Schwarzschild radii this irradiation may dominate in the lo- 
cal disc energy budget. The ratio between the irradiation 
fiux due to the warm absorber and the locally dissipated 
fiux roughly saturates at distances larger than the outer ra- 
dius of the warm absorber, and this ratio depends on the 
global warm absorber parameters Ttot, dmax, and /3. 

Since the effect of self-irradiation becomes again im- 
portant in the outermost part of the disc (a few thousands 
Schwarzschild radii) the relative importance of the two ef- 
fects depends on the adopted model parameters. 

3.3 IR/optical/UV spectra of irradiated accretion 
discs 

Examples of the irradiated accretion disc spectra are shown 
in Fig. |S| The reference model, without irradiation, is 
marked with dotted line. The spectrum is rather hard (blue) 
in the optical band. The self-irradiation of the disc changes 
the spectral slope only in the IR band, at wavelengths above 
3 /im. 

The indirect irradiation, however, may change the opti- 
cal spectrum considerably. We see that assuming very slow 
decrease of the warm absorber density with distance we ob- 
tain spectra which are practically flat in the optical band. 
The disc spectrum, however, hardens again in the UV band 
so the effect of irradiation looks like separate spectral com- 
ponent. Models with larger /3 still give the spectra consid- 
erably flatter in the optical range than the standard disc 
but the overall spectrum is practically smooth. The spectral 
range of this flattening depends on the adopted extension of 
the warm absorber. 

The amplitude of the effect depends essentially on the 
warm absorber optical depth. If the warm absorber depth 
Ttot is of the order of 0.1 (corresponding to the hydrogen 
column Nh ~ 1.5 x 10^^ cm~^) or less, no effect is visible 
and the intrinsic disc spectrum is unmodifled. 

For the model of the warm absorber adopted in Fig. |S1 
the spectral slope in the optical band ao is equal to -0.88, 
-0.56 and -0.20 for /3 = 0, 0.75 and 1.5. Here the convention 
Fi. oc is used. It compares interestingly with the ob- 
served slope: Laor et al. (1997) give the value Uo = —0.36, 
with dispersion of 0.22, at the basis of the analysis of several 
PG quasars, and ao = —0.44 was given by Vanden Berk et 
al. (1991) for quasars from SDSS survey. 

The assumed extension of the warm absorber also plays 
an important role. The dependence is the strongest if the 
gas density is uniform (middle panel in Fig. |HJ . If the outer 
radius is small only the UV part of the spectrum is modified 
while for an extended medium the modification of the UV 




Figure 6. The dependence of the irradiated accretion disc spec- 
trum on the warm absorber parameters: f3 (upper panel), Rmin 
and Rmax (middle panel) and 6 (lower panel). Parameters of the 
reference model: m = 0.75, M = 3 X W'^Mq, not = 0.6, /3 = 0, 

Rmin = WORschm, Rmax = 5000Rschm, dmax = 7r/3. Dotted 

line shows the disc spectrum without irradiation, and the long- 
dash-dot line shows the effect of the direct irradiation alone. 



part is weak while spectra in the optical/IR band are now 
much flatter. The affected spectral region is more narrow if 
Rmax /Rmin IS close to 1 and broader if this ratio is large. 

The opening angle of the cone occupied by the warm 
absorber (lower panel in Fig. [SJ changes mostly the overall 
normalization of the incident flux, similarly to the total op- 
tical depth. If this opening angle is too small (Omax < 10°) 
the presence of the warm absorber has no practical effect on 
the disc spectrum. 
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We also considered other types of the warm absorber 
density distributions. 

In one of these special cases we assumed that the density 
scales with the radius as p oc sin^[{R/ Rmin)^'^\- Such a 
distribution of density waves imitates episodes of enhanced 
and reduced outflow. However, the resulting spectra were in 
practice very similar to those obtained from our power law 
models for /3 = 0. Integration over volume effectively smears 
out any distribution inhomogeneities. 

In another of these special cases we assumed the geome- 
try of the warm absorber in the form of the empty cone cen- 
tered at the opening angle of 45°, as advocated for quasars 
by Elvis (2000) . In this case the fraction of the outgoing pho- 
tons intercepted by the warm absorber is very small (model 
described as 'empty cone' in the lowest panel in Fig. |SJ| and 
the change in the spectrum relatively small. 



4 EXAMPLES OF SPECIFIC SOURCES 

We choose two interesting sources for detailed modelling. 
Both are high Eddington luminosity objects and belong 
to the Narrow Line Seyfert 1 class. We do not expect 
an extended inner optically thin flow in this type of ob- 
jects, which are analogs of the galactic X-ray sources in 
their soft state (Pounds et al. 1995). Broad band spec- 
tra of both sources were extensively studied. One of those, 
RE J1034-I-396 (Puchnarewicz et al. 1995, Puchnarewicz et 
al. 2001), was already considered as the best case of irradi- 
ated accretion disc (Soria & Puchnarewicz 2002). The other 
one, PG 1211-1-143, is argued to show the presence of the 
warm absorber of considerable optical depth ( Pounds et al. 
2003, Gierlihski & Done 2004b). 

4.1 RE J1034+396 

The observational data for RE J1034+396 (z=0. 04214) was 
shown by Soria & Puchnarewicz (2002) and it was kindly 
provided by Liz Puchnarewicz. The nature of the optical/UV 
component is not clear; the spectra do not seem to show the 
narrow absorption lines typical for starlight (Puchnarewicz 
et al. 1995). The source is surprisingly stable in comparison 
to many other NLSl galaxies (Puchnarewicz et al. 1998) 
although hints of variability (in timescales of a fraction of 
a day) are seen in the EUVE lightcurve (Halpern, Leighly 
& Marshall 2003). Warm absorber features were not dis- 
covered so far in this source. However, the source was ob- 
served by Chandra only for 5 ksec, and by XMM only for 
15 ksec (according to http://heasarc.gsfc.nasa.gov/ I which 
is not enough to detect a highly ionized warm absorber. For 
example, no warm absorber was seen in ASCA data and 
XMM for Ton S180 and even during the preliminary anal- 
ysis of the Chandra data for the same source but detailed 
analysis of the 80 ksec Chandra data showed the presence 
of narrow absorption features (see Rozahska et al. 2004). 
Therefore, the existence of highly ionized warm absorber in 
RE J1034-h396 is not excluded. 

Our model of an irradiated accretion disc well represents 
both optical/UV continuum and the soft X-ray spectrum 
(Fig. 0. The mass of the black hole, M = 6.3 x 10^ 
is consistent with the limits 0.6 — 3.0 x 10^ Mq obtained 
by Soria & Puchnarewicz (2002). The object is accreting 



44 




log y[Hz] 



Figure 7. The observational data points for RE J1034-I-396 and 
the irradiated disc model (continuous line). Parameters of the 
model: rn = 0.84, M = 6.3 X lO^M©, not = 0.6, (3 = 0.0, 

Rmin = ^^Rschwj Rmax = 7000Rschm) Smax = 7r/3. Dotted 

line shows the disc spectrum without irradiation. Adopted cos- 
mological model: Ho = 75 km s^^Mpc"'^, qo = 1/2. 



still below the Eddington ratio but close to it (m = 0.84). 
The warm absorber responsible for scattering of the inner 
disc emission must be quite extended {Rmax = 7000Rschw)- 
Such a scattering region would smear most of the intrinsic 
flux variability happening in the timescale shorter than a few 
hours. Longer trends, if present in the intrinsic flux, should 
be also seen as corresponding variations in the optical/UV 
spectrum. 

ASCA and Beppo-SAX data did not show any clear ef- 
fect of the warm absorber (Soria & Puchnarewicz 2002) so 
the warm absorber in this source, if indeed present, must be 
completely ionized. However, the presence of strong outflow 
is rather expected for sources so close to the Eddington ratio. 
The soft X-ray spectrum of this objects is well fitted if the 
disc emission is comptonized by a medium at the tempera- 
ture r ~ 15 keV and optical depth r = 1.3 under assump- 
tion of the plane-parallel geometry (Soria & Puchnarewicz 
2002). This Comptonizing medium may be identified with 
such outflow, as suggested by Pounds et al. (2003) in the 
context of PG 1211-f 143. 

Very low value of the parameter l3 of the warm ab- 
sorber results from the spectrum being very flat in i/F^ plot. 
If a significant internal reddening is present in the source, 
the actual slope may be positive, and /3 larger. Mason et 
al. (1996) suggested the possible reddening Eb-v ~ 0.5 in 
RE J1034-f 396 at the basis of the Ha to Hp ratio of 5, much 
higher than expected for case B recombination. 

The small blue bump (Balmer continuum-l-blended iron 
lines) does not show up clearly in RE J1034-I-396. 
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Figure 8. The observational data points for PG 1211+143 and 
the irradiated disc model (continuous line). Parameters of the 
model: m = 0.75, M = 3.0 X IO^Mq, not = 0.6, f3 = 0.75, 

Rmin = WORschm, Rmax = 5000i?Schtu: dmax = 7r/3. Dotted 

line shows the disc spectrum without irradiation. Adopted cos- 
mological model: Ho = 75 km s^^Mpc"'^, qo = 1/2. 



4.2 PG 1211 + 143 

The observational data for PG 1211+143 (z=0.08090) were 
taken from Elvis et al. (1994). This source is strongly vari- 
able and the presented data are representative for the pe- 
riods when the source is bright (see e.g. Janiuk, Czerny & 
Madejski 2001). 

The shape of the UV continuum in this object clearly 
indicates the presence of the small blue bump due to the 
Balmer continuum and blended iron lines. Such a compo- 
nent is present in most spectra of Seyfert galaxies. There- 
fore, we modelled this component and added to the resulting 
disc spectrum. We adopted the spectral shape of the Balmer 
continuum and the shape of the blended iron line component 
after Wills, Netzer & Wills (1985), and we assumed equal 
amplitude of the two components (Neugebauer et al. 1987) . 
The overall amplitude of this component was assumed in 
such way that the component constitutes 37 % of the total 
spectrum at 3000 A and 9 % at 1890 A. The required nor- 
malization is rather typical for AGN. For example, in NGC 
5548 the small blue bump accounted for up to 10% of the 
flux at 1840 A and 30-40% of the flux at 2670 A (Peterson 
et al. 1991, Maoz et al. 1993). The resuh is shown in Fig.|Hl 

The mass of the black hole in our model [M — 3.0 x 
lO^M©) is consistent with the limits obtained by Kaspi et 
al. (2000) from reverberation {Mmean = 4.05l'J |^ x IO'^Mq, 
Mrms = 2.36lo 70 ^ Mq) , as well as with the resuhs 
based on the X-ray variability (log M = 7.0 ± 0.7 in Janiuk 
et al. 2001, M = 8.1 x IO^Mq in Nikolajuk, Papadakis & 
Czerny 2004). The source is accreting close to the Edding- 
ton ratio (m = 0.75). The parameters of the warm absorber 



are consistent with requirements obtained by Pounds et al. 
(2003). Gierlihski & Done (2004b) presented a good rep- 
resentation to the broad band data, including X-ray data 
from XMM-Newton, by postulating an outflowing strongly 
ionized medium with velocity dispersion v/c ~ 0.2, which 
Comptonized the disc spectrum. The value of the black hole 
mass used in this analysis (M = 6.5 x lO^M©, from Boro- 
son 2002) was somewhat higher than in our model and the 
Eddington ratio correspondingly higher (1.14 from Boroson 
2002) since disc irradiation was not included in their analy- 
sis. 

The extension of the warm absorber should greatly re- 
duce any optical variability in this source at timescales below 
~ 5 years. 

4.3 Blue radio quiet quasars 

Large sample of quasars obtained in course of completing 
Sloan Digital Sky Survey (hereafter SDSS) allowed to study 
the optical/UV slope distribution. The sources were divided 
into four classes of increasing 'reddening' in their spectra 
(Richards et al. 2003) and the composite spectra had the 
slope (measured between 1450 A and 4040 A) of -0.25, -0.41, 
-0.54 and -0.76, correspondingly. This sequence of spectra 
was well explained assuming that the effect is due to the 
increasing amount of reddening by amorphous carbon dust 
(Czerny et al. 2004). In this interpretation no irradiation of 
the disc seemed to be needed. However, the spectrum of the 
composite 1 of Richards et al. (2003), unmodified by dust, 
was not fully satisfactorily represented by a bare disc - clear 
flattening of the spectrum below log u ~ 14.9 was visible. 

We can model this flattening assuming the presence of 
the warm absorber. An example of solution is shown in 
Fig. 1^ The model shows flattening at the longest wave- 
lengths, as requested. However, we did not find a fully sat- 
isfactory representation for this high quality data - the the- 
oretical spectrum is never quite as flat as requested by the 
data of Richards et al. (2003) . 

An alternative possibility exists - Vanden Berk et al. 
(2001) attributed the change of the spectral slope in the 
quasar composite spectrum at ~ 5000 A to starlight con- 
tamination, at the basis of the detection of stellar absorption 
lines. Starburst activity may be intrinsically connected with 
an accretion process. A number of recent papers discuss the 
issue of the stars originating/existing very closely to the ac- 
tive nucleus and interacting directly with the accretion discs 
(e.g. Collin & Zahn 1999, Nayakshin 2004 and the references 
therein). Such a nuclear cluster may contribute to the overall 
spectrum. An exemplary spectrum of such a nuclear cluster 
(see Figure 1. of Nayakshin 2004) peaks at z/ « 4 x 10^"* Hz. 

We attempted to model the starlight contribution, 
taking the spectrum of the galaxy M31 as a template. 
Normalization of this contribution was a free parameter. 
Such a method was justified by the analysis performed by 
Wamsteker et al. (1990). The result is shown in Fig.|51 

The model with starlight contribution did not repro- 
duce the quasar spectrum much better than the model of 
irradiated disc. It is possible that the use of M31 as tem- 
plate is not appropriate (see e.g. Kotilainen & Ward 1994 
for the discussion of coupling between an AGN and its host 
galaxy), and better model of the circumnuclear cluster is 
needed. 
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Figure 9. The composite 1 quasar spectrum from SDSS 
(Ricliards et al. 2003; tiiin continuous line), disc spectrum with- 
out irradiation (dotted line) and disc illuminated due to the pres- 
ence of the warm absorber (thick continuous line). Parameters 
of the model; m = 0.15, M = IO^Mq, T„,ax = 0.6, = 0.75, 
Rmin = 5Q0Rschw, Rmax = lOORschw, Omax = 7r/3. Dashed 
Une shows the non-irradiated disc with starlight. 



If indeed all flattening at long wavelengths is due to 
starlight it means that the outflow of highly ionized material 
is either inefficient (low optical depth) or coUimated. 



5 DISCUSSION 

The emission from a standard accretion disc was found to 
be a very natural model of the optical/UV spectra of AGN 
(Shields 1978, Czerny & Elvis 1987, Wandel & Petrosian 
1988, Sun & Malkan 1989). However, spectral fits usually 
requested the presence of an additional underlying power 
law component of unspecified nature since the spectra were 
flatter than predicted by models, and this power law did not 
even seem to be connected to the power law component seen 
in X-ray band (Thompson 1995). 

This serious weakness of the disc model caused interests 
in other models of the optical/UV bump like optically thin 
emission (Barvainis 1993). However, detection of a Balmer 
edge in a polarized spectrum of a quasar Ton 202 by Kishi- 
moto, Antonucci & Blaes (2003) provided a direct evidence 
in the favour of the optical/UV emission originating in an 
optically thick medium (a disc). 

However, the discrepancy between the observed opti- 
cal/UV spectrum and the predictions by the canonical mod- 
els remains. It is not surprising since several effects should 
(in principle) be included in modelling optical/UV spectra 
of AGN: 

• extinction due to circumnuclear dust 



• cold disc evaporation in the innermost part and the 
transition to a hot flow 

• contribution of starlight 

• contribution of small blue bump (Balmer continuum 
and blended iron lines) 

• irradiation of the cold disc by the X-ray emission from 
the hot optically thin plasma 

• Cornptonization of the disc emission by a hot plasma 

• irradiation of the disc by scattered disc emission 

• self-irradiation of the cold disc 

• proper local models of irradiated disc atmosphere 

• non-stationarity of accretion flow and of X-ray genera- 
tion. 

Complexity of the problem is depressing, and in many cases 
it is difficult to disentangle various effects. However, some 
help comes from analyzing broad band spectra since infor- 
mation about the X-ray emission provides constraints to fac- 
tors modifying optical/UV band. 

In the present paper we model the scattering of the 
disc radiation by an extended ionized medium and redirect- 
ing its emission back again towards the disc. The observa- 
tional evidence of such a medium comes from the studies 
of warm absorbers in Seyfert galaxies (e.g. Ashton et al. 
2004 and the references therein) and Broad Absorption Line 
(BAL) quasars (e.g. Gallagher et al. 2004 and the references 
therein). Observations of these objects give direct measure- 
ments of the column density of the partially ionized material, 
being less sensitive to the fully ionized plasma. Polarimetry 
observations of Seyfert 2 galaxies with hidden Broad Line 
Regions (Antonucci & MiUer 1985; Tran 2001) and of BAL 
QSO (e.g. Lammy & Hutsemekers 2004) give in turn direct 
information about the 'scatterers' but at unspecified (rather 
large) distances from the central black hole. Estimates of the 
total column density of the material along the line of sight 
obtained from the UV studies are frequently inconsistent, 
and the amount of fully ionized material can be frequently 
underestimated. The detection rate of the warm absorber 
decreases with an increase of the object accretion rate (e.g. 
Laor et al. 1997), BAL features are seen in bright QSO. 
Although the connection between BAL outflow and warm 
absorbers is not clear they may represent basically the same 
outflow but of different ionization state and/or seen at dif- 
ferent inclination angles. 

This material can be traced, however, by modelling the 
optical/UV spectra. We show that if the warm absorber has 
considerable optical depth Ttot > 0.1 (or Nh > 10^^ cm~^), 
and if this medium is located relatively close to the central 
region we can expect significant change in the optical/UV 
spectrum of the accretion disc in such source due to the 
scattering of the fraction of radiation by this warm absorber 
and redirecting a fraction of this emission back towards the 
disc. The effect depends not only on the total optical depth 
of the warm absorber along the line of sight but on the whole 
density distribution. 

We considered in detail the case of a disc irradiated 
due to the presence of the warm absorber but we neglected 
disc evaporation and irradiation by hard X-ray source. We 
showed that the direct self-irradiation of the disc in this 
case is not very efficient but the indirect irradiation due to 
the scattering by the warm absorber may fiatten the optical 
spectrum considerably. 
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In our opinion the model directly applies to high ac- 
cretion rate AGN, m > 0.2 like quasars or Narrow Line 
Seyfert 1 galaxies. In those sources the disc extends practi- 
cally down to the marginally stable orbit (Pounds et al. 1995, 
Sobolewska, Siemiginowska & Zycki 2004), and an extended 
ADAF does not form. The hard X-ray source contains a 
small fraction of the source bolometric luminosity, as the 
optical/UV/soft X-ray Big Blue Bump strongly dominates 
the overall spectrum. Therefore, the irradiation by hard X- 
rays is energetically inefficient although it is important from 
the point of view of the reflection component formation. 
Additionally, if the hard X-ray emission comes from coronal 
flares and it is confined to the direct vicinity of the disk, our 
assumptions of emission comming from the disk surface are 
roughly satisfied even for this component. 

We obtained satisfactory representation of the data for 
two NLSl galaxies, RE J1034-h396 and PG1211-I-143. Blue 
quasar composite 1 spectrum of Richards et al. (2003) is 
almost well represented by a standard disc, with some flat- 
tening only above A ~ 5000 A. We modelled this flattening 
also as the effect of irradiation due to warm absorber but 
the result was not quite convincing and the contribution of 
starlight seems to be most probable explanation. 

We neglected other possible effects from our list al- 
though intrinsic reddening as well as Comptonization can 
strongly bias the observed FUV spectral index (Shang et al. 
2004). 

The irradiation of the disc due to the presence of the 
warm absorber may also be important in Seyfert 1 galax- 
ies. However, typical Seyfert 1 galaxies may have their discs 
truncated at larger radius, and a geometrically thick inner 
flow develops (e.g. Poutanen, Krolik, & Ryde 1997, Chiang 
& Blaes 2003), the bolometric luminosity of the hard X- 
ray component is larger, and the source is most probably 
located at large distance from the equatorial plane, either 
corresponding to the disk truncation radius (in ADAF mod- 
els), or to the position of a base of jet (Henri & Pelletier 1991; 
Martocchia, Matt & Karas 2002). It enhances the direct self- 
irradiation by the central X-ray emitting region (Chiang & 
Blaes 2003). So it is possible that in Seyfert galaxies both 
direct X-ray irradiation and indirect irradiation due to the 
warm absorber scattering can be equally important, but the 
model would include additional parameters specifying the 
location of the X-ray emitting plasma. 
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